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(54) Magnetic resonance imaging 

(57) In magnetic resonance imaging apparatus em- 
ploying magnetic gradient fields in a phase-encode and 
in a read-out direction for spatially encoding excited MR 
active nuclei in a region of interest of a patient, in which 
a reduced number of readings in the read-out direction 
is taken, thereby creating an aliased reduced field of 



view image, at least two r.f. receive coils are used to- 
gether with sensitivity information concerning those 
coils in order to unfold the aliased image to produce a 
full image while taking advantage of the reduced time of 
collection of data. In accordance with the invention, sen- 
sitivity information is collected at a lower resolution than 
that at which the image information is collected. 
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Description 

[0001] This invention relates to magnetic resonance 
(MR) imaging. 

[0002] The invention is particularly concerned with re- 
duction in the time needed to collect data for leasing a 
region of interest of a patient. 

[0003] For example, it is sometimes desired to view 
the progress of a contrast agent in the bloodstream 
through a region, in order to help distinguish between 
benign tumour, malignant tumour and unaffected tissue, 
and rapid imaging is required to monitor the flow of con- 
trast agent to the various sites. A malignant tumour may 
require a greater supply of blood and could show up in 
this way. 

[0004] A prior art magnetic resonance imaging appa- 
ratus is shown in Figure 1 . A patient 1 (shown in section) 
is slid axially into the bore 2 of a superconducting mag- 
net 3, and the main magnetic field is set up along the 
axis of the bore, termed by convention the Z-direction. 
Magnetic field gradients are set up, for example, in the 
Z-direction, to confine the excitation of magnetic reso- 
nant (MR) active nuclei (typically hydrogen protons in 
water and fat tissue) to a particular slice in the Z-direc- 
tion e.g. that illustrated in Figure 1 and, in the horizontal 
X and the vertical Y directions as seen in Figure 1, to 
encode the resonant MR nuclei in the plane of the slice. 
An r.f transmit coil (not shown) applies an excitation 
pulse to excite the protons to resonance, and a pair of 
r.f. receive coils 4, 5 pick up relaxation signals emitted 
by the disturbed protons. 

[0005] To encode/decode received signals in the Y- 
direction, the signals are detected in the presence of a 
magnetic field gradient, termed a read-out gradient, to 
enable different positions of relaxing nuclei to corre- 
spond to different precession frequencies of those nu- 
clei about the direction of the main magnetic field due 
to the influence of the gradient. The data is digitised, 
and so for each r.f. excitation pulse, a series of digital 
data points are collected, and these are mapped into a 
spatial frequency domain known as K-space (Figure 2). 
Each r.f. pulse permits at least one column of digital data 
points to be collected. 

[0006] To encode/decode the received signals in the 
X-direction, after each r.f. pulse has been transmitted 
and before data is collected with the read-out gradient 
applied, a magnetic field gradient in the X-direction is 
turned on and off This is done for a series of magnitudes 
of magnetic field gradients in the X-direction, one r.f. 
pulse typically corresponding to a different magnitude 
of gradient in the X-direction. 

[0007] On the K-space matrix shown in Figure 2, the 
columns of data points correspond to r.f. pulses followed 
by different magnitudes of phase-encode gradients. 
[0008] The field of view imaged by the magnetic res- 
onance imaging apparatus depends on the spacing of 
the data points in the phase-encode and read-out direc- 
tions, and the resolution of the image depends on how 
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far the points extend in each direction e.g. how large the 
maximum phase-encode gradient is. 
[0009] The signals received from the r.f. receiver coils 
4, 5 are subject to a two dimensional fast Fourier Trans- 
5 form in Fourier Transform processors 6, 7, to produce 
. pixelateH imofjes which are $tore.d in insane memories 
8,9. 

[0010] A slice image is shown in Figure 3. For the pur- 
poses of explanation, the symbols of a circle, a square 
10 and a triangle have been illustrated in both the object 2 
and the image in Figure 3. Figure 3 implies thatthe spac- 
ing of data points in the phase-encode gradient direction 
is sufficient to image the whole of the slice shown in Fig- 
ure 1 . 

15 [0011] Between each r.f. pulse, there is a certain min- 
imum pulse repetition time, and the collection of data 
implied by Figures 2 and 3 may therefore take too long 
to permit flows of contrast agent to be adequately mon- 
itored. 

20 [0012] One technique used to reduce the data collec- * 
tion time is to cut out, say, half the phase-encode steps 
e.g. by keeping the same maximum phase-encode gra- 
dient but omitting every other column of data. This would 
then halve the data collection time. 

25 [0013] The spacing of the data points in the phase- 
encode direction would now have doubled, so that the 
field of view in the corresponding image domain would 
have halved. (The field of view in the read-out direction 
would remain the same because the number of data 

30 points collected during read-out would remain the 
same.) The imaged area would now cover little more 
than the square illustrated in Figure 1 . This is shown in 
Figure 4. Unfortunately, MR signals from the circle and 
the triangle would still be picked up by each receive coil 

35 4, 5, and the processing is such that these regions would 
be folded back over the square, with the circle and the 
triangle being interposed. This is a problem known as 
aliasing. Separate points P v P 2 in the object would be 
imaged at the same pixel in Figure 4. 

40 [001 4] It is to enable the aliased image to be unfolded 
that two receive coils 4, 5 have been used. (Magnetic 
Resonance in Medicine 42: 952-962 (1999) - SENSE: 
Sensitivity Encoding For Fast MRl by Klaas P Pruess- 
mann, Markus Weiger, Markus B Scheidegger and Pe- 

45 ter Boesiger). With the two receive coils and suitable 
processing, the benefit of a halved data collection time 
can be retained, but the image can be unfolded. This 
technique is referred to as SENSE (sensitivity encod- 
ing). 

so [0015] The response of coil 4 to MR signals from point 
P 1 is perhaps greater than its response to MR signals 
emanating from point P 2 , and the position is perhaps the 
other way round with regard to the coil 5. Knowledge of 
the relative sensitivities of the coils 4, 5 can be used in 

55 the processor 10 to unfold the aliased image of Figure 
4 into an image having the pattern of Figure 3, although 
with a lower signal-to-noise ratio because only half the 
data has been collected. 



EP1 102 076 A2 



2 




3 

[0016] For example, coil 4 receives MR signals from 
both points P 1 and P 2 , and the combined intensity of the 
resulting pixel stored in image memory 8 (C^ depends 
on the sensitivity of the coil 4 to point (S 1n ) and the 
sensitivity of the coil 4 to point P 2 (S 12 ). The image 5 
memory 9 also has an intensity stored (C 2 ) correspond- 
ing to points P land P 2 which depends on the sensitivity 
of coil 5 to P 1 (S 21 ) and on its sensitivity to P 2 (S 22 ). 
Thus, 

C 1 =S 11 P 1 +S 12 p 2 
C 2 =S 21 P 1 + S 22 P 2 

[0017] More generally (C) = (S).(P) where (C), (S) and 
(P) are matrices. 

[0018] Provided the sensitivity terms are known, the 
intensity correspond^ i y :o each pixel ?^ and r 2 sepa- 
rately can be calculated in processor 10, for each pair 
of points of the image which are mapped onto the same 
pixel in the aliased image. The aliased image can then 
be unfolded and viewed in full on display 11 . 
[0019] To calibrate the coils 4, 5 requires a knowledge 
of the sensitivity of each coil at each point. Of course, 
the data collected at each image memory 8, 9 maps the 
sensitivity of each coil 4, 5 onto each pixel, but in a way 
which is clearly weighted by the intensity produced by 
each point of the object. 

[0020] The paper referred to suggests collecting the 
sensitivity data using a third receive coil, in the form of 
a body coil, which has a uniform sensitivity over each 
part of the object. 

[0021] To suppress the weighting produced by the ob- 
ject being imaged, it turns out that the lines of the images 
in the image memories 8, 9 have to be divided by cor- 
responding lines i.e. for the same point in K-space in the 
read-out direction, produced by the imaging with the 
body coil. 

[0022] However, for example, in the case of an object 
with the weighting shown in Figure 5, the sensitivity pro- 
file for the coil 4 could be as shown in Figure 6 and that 
for the body coil as shown in Figure 7. Figure 8 shows 
the division of one by the other and the presence of 
noise. 

[0023] This is because the body coil picks up noise 
from the whole of the section of the body being imaged, 
and the resulting sensitivity profile produced by the di- 
vision, for each line of the image, is masked by noise. 
[0024] To overcome this , the paper referred to propos- 
es to divide the response of one coil 4 (Figure 9) by the 
response of the other coil 5 (fiyuie 1C) and to dispense 
with the body coil altogether, for each line of the image. 
This is possible because it turns out that the (S) matrix 
terms only require this ratio to suppress the weighting 
effect of the object and the correction for the relative 
sensitivities produced by the two coils. 



[0025] There is now less noise in the resultant division 
(Figure 11) but a problem arises which is also inherent 
in the ratio shown in Figure 8 i.e. that the ratio of sensi- 
tivities is only present where the object is present. Noise 
is still present, but to a lesser extent, and postprocessing 
therefore has to be carried out to reduce the nqisfi. The. 
problem is that, in so-doing, a postprocessing algorithm 
has to be employed which prevents similar sharp edges 
in the edges of the object itself being smoothed. 
[0026] The invention provides apparatus for magnetic 
resonance imaging, comprising means for creating 
magnetic field gradients in a phase-encode direction for 
spatially encoding excited MR active nuclei in a region 
of interest of a patient, means for receiving r.f. signals 
from the excited M R active nuclei, means for processing 
the spatially encoded r.f. receive signals to produce a 
spatial representation of the region of interest, the mag- 
netic field gradient means being arranged to create a 
number of phase-encode gradients corresponding to a 
smaiier iieid ot view than the region oi interest, the r.f. 
receive means including a pair of receive coils and 
processing means arranged to use the r.f. signals re- 
ceived from each coil together with sensitivity informa- 
tion for the coils to produce an unfolded version of the 
aliased data corresponding to the reduced field of view, 
wherein the sensitivity information used is arranged to 
be collected at reduced resolution compared to that at 
which the image is collected. 

[0027] The invention also provides a method of mag- 
netic resonance imaging, comprising creating magnetic 
field gradients in a phase-encode direction to spatially 
encode excited MR active nuclei in a region of interest 
of a patient, receiving r.f. signals from the excited MR 
active nuclei, processing the spatially encoded r.f. re- 
ceive signals to produce a spatial representation of the 
region of interest, the number of phase-encode gradi- 
ents being created corresponding to a smaller field of 
view than the region of interest, a pair of r.f. receive coils 
receiving the r.f. signals, and the processing using the 
r.f. signals received from each coil together with sensi- 
tivity information for the coils to produce an unfolded ver- 
sion of the aliased data corresponding to the reduced 
field of view, wherein the sensitivity information used is 
collected at reduced resolution compared to that at 
which the image is collected. 

[0028] By collecting the sensitivity information at a 
lower resolution, it become unnecessary to fit or smooth 
the sensitivity information removed of its object weight- 
. ing. 

[0029] The sensitivity information may be collected by 
ratioing each of the receive coils with a body coil, or by 
ratioing the sensitivities of the two r.f. receive coils them- 
selves. These may be of the loop and butterfly type as 
claimed in our co-pending British patent application no. 
9926923.5, but could be a pair of loops or any other 
types. Further, further pairs of loops may be provided. 
[0030] Ways of carrying out the invention will now be 
described in greater detail, by way of example, with ref- 
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erence to the accompanying drawings, in which: 

Figure 1 is a partly axial sectional, partly block dia- 
gram representation, of prior art magnetic reso- 
nance imaging apparatus; 

Figure 2 is a representation of the data in K-space 
resulting from the signal picked up by each receive 
coil of the apparatus; 

Figure 3 is a representation of the image area which 
would be built up if a full number of phase-encode 
gradients were performed; 

Figure 4 is an image which would be built up if a 
reduced number of phase-encode gradients were 
to be performed; 

Figure 5 represents the object weighting of a line in 
the imago met nCiy 5, 9, 

Figure 6 represents the object weighted sensitivity 
of coil 4 for a line in the image memory 9; 

Figure 7 represents the object weighted sensitivity 
of a body coil for an equivalent line in an image 
memory connected to the body coil; 

Figure 6 represents the ratio of the object weighted 
sensitivities shown in Figures 6 and 7; 

Figure 9 represents the object weighted sensitivity 
of coil 4 for a line in image memory 9; 

Figure 1 0 represents the object weighted sensitivity 
of coil 5 for a line in image memory 8; 

Figure 11 represents the ratio of the object weighted 
sensitivities shown in Figures 9 and 10; and 

Figure 12 is a partly axial sectional, partly block di- 
agram representation, of magnetic resonance im- 
aging apparatus according to the invention. 

[0031] In accordance with the invention (Figure 12), 
reference data is collected at reduced resolution com- 
pared to that used during imaging. The reference data 
is of course collected at full field of view. Lower maxi- 
mum values of phase-encode gradient are used. The 
effect of lower resolution in the reference data, used to 
calibrate the sensitivity of the coils 4, 5 is to reduce noise 
in the reference data. The reference data, which is 
•stored in memory 12, is used in the SENSE method to 
unfold the aliased data of Figure 4. Because the noise 
is reduced in the reference data, the signal-to-noise of 
the target unfolded SENSE data, is increased. 
[0032] Various arrangements may be used. One pair 
of coils 4, 5 may be used. A loop, butterfly combination 
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as claimed in our co-pending British patent application 
no. 9926923.5 may be used. More than one pair of coils 
may be used. In each case, the reference data may be 
produced using a body coil as in Figures 6 to 8, or a 

5 body coil may be omitted, as in Figures 9 to 11 . 

[0033] To go now into the background of the invention 
in greater detail, dynamic enhanced MR mammography 
(DMRM) is a leading imaging method for identifying 
breast cancer in patients with dense breast tissue. Dif- 

10 ferentiation of malignant from benign disease is critically 
dependent on contrast enhancement profiles; cancers 
enhance the most briskly with peak enhancements at a 
maximum of two minutes. Current screening techniques 
rely on coverage of both breasts using volume acquisi- 

15 tions following bolus injection of Gadolinium-DTPA. The 
requirements for large field of view (FoV) at millimetre 
resolution with adequate temporal resolution are chal- 
lenging. 

[0034] Sensitivity encoding SENSE requires an array 

20 coii with uoii elements that hate spatially tfioii.^t Scnoi- 
tivity profiles over the region of interest. A 4 channel re- 
ceive only array coil with two pairs of nominally orthog- 
onal coils is used; each pair providing approximately 
uniform sensitivity over one breast. Only one pair of coils 

25 is shown in Figure 12. The dynamic examination con- 
sisted of full resolution volume scans with the phase en- 
coded, left-right FoV reduced by two as explained in our 
co-pending British patent application no. 9926923.5. 
The resulting degeneracy is unfolded using sensitivity 

30 data from the coils. 

[0035] To avoid misalignment between reference and 
target data and changes in coil loading, sensitivity data 
is acquired with the subject in situ. This data is contam- 
inated by a dependence on tissue signals, which can be 

35 removed by normalising using anatomical images that 
are assumed to be free of sensitivity variations. Alterna- 
tively, the approach used here uses the relative sensi- 
tivities of the coils in the array. In either case the sensi- 
tivity data requires a high signal-to-noise ratio (SNR) to 

*o avoid compromising the SNR of the final data. Filtering 
the sensitivity data can increase SNR. Artefacts in the 
sensitivity measurements from cardiac or respiratory 
motion are a critical issue in the current application. 
Non-localised artefacts tend to mask the spatial proper- 

45 ties of the coils. 

[0036] To provide high SNR reference data that is ro- 
bust against physiological motion we have reduced the 
resolution of the full FoV reference acquisitions, but 
maintained approximately constant imaging time by 

50 multiple averaging. This reverses the order of the nor- 
malisation and filtering steps, leaving the coil profile in- 
tact provided the reference data has sufficient spatial 
frequency content. 

[0037] In one embodiment, all data was acquired on 
55 a 0.5T. Picker Apollo scanner (Cleveland, OH) using a 
3D gradient echo sequence (TR/TE 10/3.75) to image 
a region 30x3ox18cm 3 . Coil sensitivity information was 
obtained with a matrix of 64x32x50, 1 0NEX and was ze- 
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ro filled to 256x256x50 prior to reconstruction. Dynamic 
data was acquired in 10 frames, each with harf left-right 
FoV and a rectangular matrix of 256x128x50, 1NEX 
(40s per frame). The data were reconstructed and un- 
folded onto a full FoV, 256x256x50 matrix. Unfolding 
was performed pair wise for the two symmetrical left- 
right coil combinations. The final images were then add- 
ed in quadrature. Image manipulation employed IDL 
(Research Systems, Colorado). 
[0038] The method was tested on phantoms, two nor- 
mal volunteers without contrast, and one normal volun- 
teer with contrast. Reference data acquired at full reso- 
lution with a single average and reduced resolution with 
multiple averages were collected and compared for 
SENSE processing. Reduced FoV SENSE processed 
data and full FoV standard images were obtained using 
the same acquisition. Images were assessed for ana- 
tomical detail, artefacts and sensitivity to contrast up- 
take. 

[GG39] Use of idii resolution single average reference 
data resulted in errors in the unfolding associated with 
cardiac and respiratory motion. The use of low resolu- 
tion, multiply averaged reference data eliminated this 
problem to produce artefact free unfolded images with 
no detectable loss of resolution. 
[0040] In this embodiment, an application of SENSE 
increases the temporal resolution of DMRM by a factor 
of two. The method has been shown to have sufficient 
sensitivity and immunity from artefacts for detecting en- 
hancement in normal tissue with a temporal resolution 
well suited for pathological uptake 
[0041] While the magnet has been described as a su- 
perconducting magnet, it could be a resistive magnet, 
and need not be of the annular configuration but could 
be an open C or H configuration, and indeed, permanent 
magnets could be used. 



Claims 

1. Apparatus for magnetic resonance imaging, com- 
prising means for creating magnetic field gradients 
in a phase-encode direction for spatially encoding 
excited MR active nuclei in a region of interest of a 
patient, means (4, 5) for receiving r.f. signals from 
the excited MR active nuclei, means (6-9, 10, 12) 
for processing the spatially encoded r.f. receive sig- 
nals to produce a spatial representation of the re- 
gion of interest, the magnetic field gradient means 
being arranged to create a number of phase-en- 
code gradients corresponding to a smaller field of 
view than the region of interest, the r.f. receive 
means (4, 5) including a pair of r.f. receive coils and 
the processing means (6-9, 10, 12) being arranged 
to use the r.f. signals received from each coil togeth- 
er with sensitivity information (12) for the coils to 
produce an unfolded version of the aliased data cor- 
responding to the reduced field of view, wherein the 



sensitivity information used is arranged to be col- 
lected at reduced resolution compared to that at 
which the image is collected. 

5 2. Apparatus as claimed in Claim 1 , in which the sen- 
sitivity information is arranged to be produced using 
a lower value of maximum phase-encode gradient 
than that used for collecting the image information. 

10 3. Apparatus as claimed in Claim 1 or Claim 2, in which 
the sensitivity information is arranged to be collect- 
ed by ratioing the sensitivities of two r.f. receive 
coils. 

15 4. Apparatus as claimed in Claim 1 or Claim 2, in which 
the sensitivity information is arranged to be collect- 
ed by ratioing the sensitivity of one r.f. receive coil 
with a body coil. 

20 5. Apparatus c.3 clalrr.cJ ^.y one of CIsiinD i to 4, 
in which the sensitivity information is arranged to be 
collected a plurality of times and averaged. 

6. A method of magnetic resonance imaging, compris- 
es ing creating magnetic field gradients in a phase-en- 
code direction to spatially encode excited MR active 
nuclei in a region of interest of a patient, receiving 
r.f. signals from the excited MR active nuclei, 
processing the spatially encoded r.f. receive signals 

30 to produce a spatial representation of the region of 
interest, the number of phase-encode gradients be- 
ing created corresponding to a smaller field of view 
than the region of interest, a pair of r.f. receive coils 
receiving the r.f. signals, and the processing using 
35 the r.f. signals received from each coil together with 
sensitivity information for the coils to produce an un- 
folded version of the aliased data corresponding to 
the reduced field of view, wherein the sensitivity in- 
formation used is collected at reduced resolution 
40 compared to that at which the image is collected. 

7. A method as claimed in Claim 6, in which the sen- 
sitivity information is produced using a lower value 
of maximum phase-encode gradient than that used 

45 for collecting the image information. 

8. A method as claimed in Claim 6 or Claim 7, in which 
the sensitivity information is collected by ratioing the 
sensitivities of two r.f. receive coils. 

50 

9. A method as claimed in Claim 6 or Claim 7, in which 
the sensitivity information is collected by ratioing the 
sensitivity of one r.f. receive c-i! v;:th p Hody coil. 

55 10. A method as claimed in any one of Claims 6 to 9, in 
which the sensitivity information is collected a plu- 
rality of times and averaged. 
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